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Summary. — We provide a brief summary of recent developments in QCD cal- 
culations in and beyond fixed-order perturbation theory for observables in e'^e~ 
annihilation as well as hadron collisions. 



1. Introduction 

The interplay between physics at hadron colliders and that at e^e^ machines has 
traditionally been of great significance in furthering our understanding of elementary 
particles and their interactions. One can for example point to the specific case of the 
discovery of the Z boson at a hadron collider [1] which was then followed by high precision 
phenomenology at LEP which helped to establish firmly the standard model of particle 
physics, the current theory of elementary particles beyond which any discoveries are still 
to be made. 

This tradition is set to continue with the strong expectation that the LHC will lead to 
the discovery of the Higgs boson or help to clarify the Higgs sector as well as enabling the 
discovery of physics beyond the standard model. The extremely high energy hadronic 
collisions at the LHC, which make it a powerful discovery machine, however come with 
a price which takes the form of a more complicated initial state (protons rather than 
elementary particles) and complications concerning non-perturbative effects such as beam 
remnant interactions (the underlying event) and pile-up which threaten to limit the 
theoretical precision that one may be able to obtain. The most precise determination of 
the parameters of the new physics such as masses and couplings would probably require 
a cleaner environment such as a high energy e+e" future linear collider. 

Nevertheless as the Tevatron experience has to an extent confirmed, calculations in 
perturbative QCD will have a strong role to play in the physics program of the LHC, 
particularly with regards to estimating accurately backgrounds to new physics. To this 
end significant effort has been devoted in the past years to develop QCD calculations 
specifically for important LHC processes in the discovery context. Moreover given the 
vast scale hierarchy inherent in high energy hadron collider physics (with scales ranging 
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from the TeV range centre-of-mass energy y^, through typical jet transverse momenta 
Pt, the masses of electroweak scale particles down to the few GeV scales associated to 
non-perturbative physics) it is clear that techniques involving summation of large loga- 
rithms in scale ratios would also be important in maximising the theoretical accuracy one 
may be able to achieve. The introduction of new and faster infrared and coUinear (IRC) 
safe jet algorithms and a systematic understanding of perturbative and non-perturbative 
properties of jets and jet substructure is also a rapidly developing and vital part of the 
current and future LHC physics program. 

At the same time, as should be clear from the preceding discussion, furthering the 
precision of QCD calculations for e^e~ annihilation remains of continued importance 
for future phenomenology as well as remaining a simpler learning and testing ground 
for QCD practitioners. In this context the development of next-to-next to leading or- 
der (NNLO) predictions and taking resummed computations from the state of the art 
next-to-leading logarithmic (NLL) level through to NNLL accuracy as well as possi- 
bly improving the current theoretical understanding of non-perturbative effects such as 
hadronisation corrections will all play an important role. 

In what follows below we present a brief summary of what we perceive to be some of 
the main developments and recent progress in QCD calculations for both hadron colliders 
and e^e~ machines. It is impossible due to page limitations to adequately cover all the 
relevant progress that has been made in the past few years and thus the selection of 
topics/references below is far from complete. We shall aim to discuss briefly the progress 
in fixed-order perturbative computations as well as all-order resummations both in the 
hadron collider and the e+e^ context, mention the status of Os measurements and discuss 
progress in the definition and understanding of jets and their properties in and beyond 
QCD perturbation theory. 

2. QCD at fixed order 

Observables that have the property of infrared and coUinear (IRC) safety can be 
calculated as an expansion in the strong coupling as using perturbative techniques based 
on the evaluation of Feynman graphs. By an IRC safe observable one essentially means 
the following: Let ©„ = 0(^1,^2, • • 'Pn) denote the value of the observable O due to a 
configuration involving n partons with momenta pi,P2' • • Pn- Now consider adding an 
extra parton with momentum Pn+i- In the soft limit that the energy En+i — > (with 
-Bi, • • • , En held finite) or the limit that Pn+i — > Pi where i — 1,- ■ ■ ,n i.c the limit in 
which the three-momentum Pn+i is parallel to any of the three-momenta pi (coUinear 
limit) IRC safety implies that independent of n, On+i — > On- IRC safety ensures that 
real-virtual cancellation of divergences occurs and hence that finite results are obtained 
in perturbation theory. 

For a simple observable of the above kind, V, involving a single hard scale Q^, we can 
then write the perturbation expansion as 

(1) v^Y.Cn(%)a:{^,% 

where C„ are perturbatively calculable coefiicients, Q is the hard scale of the process 
and /i an arbitrary renormalisation scale, which however should be chosen to be of order 
Q to avoid large logarithms in Q^ / fp. The dependence on /i would in fact cancel if one 
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Fig. 1. - An illustration of the scale uncertainty reduction with the order of the perturbative 
estimate for the case of the rapidity {Y) distribution for inclusive Z production at the LHC. 
Figure taken from Ref. [3]. 



were able to compute the observable to all orders exactly but in practice one is able to 
evaluate only a few terms in the above sum. The residual ^ dependence in a calculation 
truncated at v}^ order in as is of the order of uncalculated O (a""*"^) terms. Thus scale 
dependence is usually taken as a measure of the influence of uncalculated higher orders 
and hence the theoretical accuracy of a given prediction. {}) 

Generally speaking leading order (LO) calculations are too crude to be considered 
reliable estimates for most collider observables. NLO calculations on the other hand 
may be expected to be correct, broadly speaking, to within order 10 percent while NNLO 
calculations represent high precision and as a rule of thumb ought to be accurate to within 
a few percent or so. (^) An illustration of this is provided in Fig. 1 where one notes the 
progressive reduction in scale uncertainty with the increasing order of the perturbative 
evaluation for the case of the inclusive Z rapidity distribution for the LHC. 

For reliable estimates of backgrounds to LHC processes with new physics it would 
thus appear that at least NLO accuracy is a must. For up to the production of three 
jets at hadron colliders NLO calculations encoded in the program NLOjet++ have been 
available for some time [4] . However many of the relevant discovery processes involve high 
multiplicity final states with similar backgrounds involving e.g multiple hard final state 
jets for which it is much less straightforward to obtain NLO estimates. At present the 
current state of the art for NLO computations at hadron colliders is for 2 -> 4 processes 
such as a ttbb final state relevant for Higgs production and decay in association with a tt 



(}) One should be aware that the scale dependence may in cases not be a reliable estimate of 
the true size of higher orders. For example if new hard scattering channels open up in higher 
perturbative orders varying scales in a lower order contribution cannot be expected to estimate 
the size of such new contributions. 

(^) There are exceptions to these broad statements which for instance only apply to observables 
not afflicted by multiple disparate hard scales. For explicit counter examples see for instance 
Ref. [2]. 
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pair [5, 6]. Similarly NLO calculations to W + 3j [7, 8] and Z'^ + 3j [9] have been recently 
computed. A significant development in the computation of NLO corrections has been the 
advent of unitarity based calculational methods alongside traditional Feynman-diagram 
techniques. A pedagogical review and further references can be found in Ref. [10]. The 
automation of NLO computations is also an important step towards the calculation of 
several different collider processes. The automation of both real radiation terms [11, 12, 
13] and virtual corrections [14, 15, 16, 17, 18] has been achieved in the past few years, 
for NLO corrections. 

As far as NNLO calculations are concerned only a few processes are known to such 
accuracy. For instance for hadron coUisons fully exclusive NNLO corrections to vector 
boson production have been computed [19, 20] while for the case of e+e^ annihilation 
similar calculations have been performed using the method of antenna subtraction for 
the case of e"'"e~ — > 3j which has enabled a more accurate determination of as from data 
on LEP event shape variables [21, 22, 23]. 

Having briefly summarised the state of the art for QCD calculations at fixed-order we 
shall turn our attention to those observables where the involvement of more than one per- 
turbative scale forces us to go beyond fixed-order perturbation theory using resummation 
methods. 

3. QCD beyond fixed-order 

As mentioned above there exist several observables of phenomenological interest where 
multiple scales (typically the process hard scale and other scales introduced due to ob- 
servable definition) play an important role. For such observables, the classic examples 
of which remain event or jet shape variable distributions [24], one has to consider the 
role of large logarithms in scale ratios and examine the possibility to resum these to all 
orders at a given logarithmic accuracy. 

To be more explicit consider the distribution in some shape variable t in say e+e^ 
annihilation: 

(2) i^^y-ic,Jln2"-ii + ... 

n 

The above behavior reflects the double-logarithmic enhancement of the shape cross- 
section due to soft and collinear gluon emissions while the ellipsis denote less singular 
terms some of which also need to be accounted for for phenomenological purposes. This 
result is clearly divergent and unphysical at small t which reflects the inadequacy of 
fixed-order predictions in that region and hence the need for resummation. 

On resummation, for those variables that have the property of exponentiation [25] 
one can write a result of the form 

, s 1 d(T d n r. lv,2 1 

(3) ^e^Cp^sin -_!_... 

a dr dr 

which generalises with account of running coupling and less singular terms into the form 
(L = lnl/r): 

(4) --T ^ ~r cxp [Lgi{asL) + g2{asL) + asgsiasL) H ] . 

a dT dT 
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Fig. 2. - Figure illustrating the comparison between various levels of resummation and results 
from PYTHIA for the thrust distribution in e^e~ annihilation for Q = 91.2 GeV (left) and 
Q = 1 TeV (right). Data from LEP are also shown in the former case. Figure taken from 
Ref. [26] 



In the above result the leading and next-to-leading logarithmic (NLL) terms are rep- 
resented by the functions gi and g2. The current state-of-the art for most observables 
at any collider process is NLL accuracy in the resummed exponent. The NNLL func- 
tion ga is known only for some select variables amongst which are the thrust and heavy 
jet-mass distribution in e^e^ annihilation (in fact computed in the framework of soft- 
coUinear effective theory to N'^LL accuracy [26]) and for hadron collisions the Drell-Yan 
and Higgs transverse momentum (Qt) distribution (see for instance Ref. [27] and ref- 
erences therein). Most recently for the Drell-Yan case results have also been obtained 
including NNLL accuracy for the new ax and 4>* variables measured by the DO collab- 
oration [28] which broadly speaking are in good agreement with the data even without 
inclusion of non-perturbative effects [29, 30]. 

For e+e^ annihilation the role of resummation in ensuring precision phenomenology 
has been clear for a long time [25] . Consider as a recent example the comparison between 
various levels of resummation, event generators and e+e" event shape data depicted in 
Fig. 2. At the Z peak it appears that there is excellent agreement between PYTHIA 
(at hadron level) and data. Moreover the PYTHIA (parton level) result appears rather 
closer to the N'^LL (4"* order) curve than to the LL result which is where one may expect 
it to be. That this is an effect which arises due to uncontrolled sub-leading terms and 
the tuning procedure inherent in PYTHIA is revealed by going to Q = 1 TeV, where 
for example subleading effects would be inconsequential, PYTHIA is much closer to the 
LL rather than the N'^LL result. It has hence been observed in Ref. [26] that using 
LL MC generators may potentially lead to a significant underestimate of certain QCD 
backgrounds at a future ILC (at about the few tens of percent level) . 

While accurate resummed predictions have been an important requirement in say the 
determination of a^ from LEP event shapes, they are also in principle of great value for 
jet production in hadron collisions in terms of improving perturbative accuracy. However 
the more complex hadronic environment at a hadron collider makes all-order resumma- 
tion a rather delicate affair. For instance care has to be taken in constructing observables 
such as event shapes to avoid contamination from beam remnants by constructing suit- 
ably central event shapes which then have the property of being non-global [32, 33]. 
Since the non-global single logarithms cannot be computed beyond the large Nc limit, 
in order to ensure full NLL accuracy for observables such as event shapes in hadronic 
dijet production, the observables have to be further modified in such a way so as to 
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Fig. 3. - Resummed predictions for global hadron collider event shapes compared to results from 
HERWIG. Figure taken from Ref. [31] 



ensure globalness, such as those variables studied in Ref. [31]. A yet more troublesome 
issue is the contamination as a result of effects such as pile-up which can potentially 
override the eventual accuracy vifhich can be achieved via theoretical methods such as 
resummation. It is thus desirable to seek variables that are less prone to such effects in 
order to test resummed calculations hadron collider observables. Predictions for several 
hadron collider event shape variables as reported in Ref. [31] are shown in Fig. 3. In some 
cases some discrepancy with corresponding results from leading-log and leading colour 
event generators such as HERWIG can also be noted. For more detailed comments on 
the role of tuning and the shower parameters in such comparisons we refer the reader to 
the comments in Ref. [31]. Detailed phenomenological studies for hadron collider event 
shape variables are currently in progress [31, 34]. 

Aside from a limited number of global event shapes and observables such as suitably 
defined dijet azimuthal correlations [35] as well as Drell-Yan Qt spectra, one may try 
to study via resummation other observables involving for instance jet-definition and the 
application of a jet algorithm. As an example of this one can point to the case of jet 
masses and shapes for high pr jets at the LHC which are relevant in identifying the 
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Fig. 4. - Figure displaying the as values extracted from various QCD studies alongwith the world 
average value (dashed vertical line) and error (yellow band) . Taken from the comprehensive 2009 
review [37]. 



origin of a jet as being initiated by a QCD process (quark or gluon jet) or say by the 
decay of a boosted heavy particle. The QCD jet mass distribution for example would 
receive logarithmic enhancements ~ a^ In -^p- where Pt is the transverse momentum 
and Mj the jet mass, with R being the jet radius. Since at the LHC we will encounter 
jets with Pf in the TeV range, the role of such logarithmic terms can be expected to 
be substantial even up to jet masses near the electroweak scale. The resummation of 
such logarithms while being immensely desirable from the standpoint of perturbative 
accuracy however has complex issues mainly to do with the role of non-global logarithms 
and jet algorithms and was recently discussed in Ref. [36]. While a very high formal level 
of precision in such cases is essentially ruled out it should still be possible to develop 
resummation formulae that capture the numerically dominant terms in the result to 
sufficient accuracy for phenomenological purposes. 

We conclude this section by presenting in Fig. 4, the current status of ag determina- 
tions taken from Ref. [37]. The 2009 value for the world average for as{Mz) was reported 
as 0.1184 ±0.0007. The individual contributions from various QCD observables used for 
as extraction are also shown in Fig. 4. 

4. Progress in jet definiton and understanding jet properties 

Although the precise definition of QCD jets may appear a detail not necessarily di- 
rectly relevant to progress of high order QCD calculations discussed in the major part 
of this review, it is in fact the case that such calculational developments need to be 
supported by suitable IRC safe jet definitions. In other words higher order perturbative 
estimates for jet cross-sections and differential distributions only make sense when an 
infrared and coUinear safe jet algorithm is used in jet definition. Although in many cases 
of interest such as inclusive jet pr spectra the IRC unsafety of a given jet algorithms 
may only appear at a relatively high order, for several LHC processes involving large 
multiplicity of final state jets (say as backgrounds to a new physics process) the IRC 
unsafety may appear already at leading order invalidating any level of perturbative ac- 
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curacy [38] . Likewise it is not meaningful to compute all-order resummcd predictions for 
quantities that will diverge at any fixed order due to the algorithm in use. This require- 
ment coupled with experimental and practical considerations (speed of the algorithm 
for high multiplicity hadronic final states) make the definition of jets a non-trivial task. 
Fortunately there now exist several different practically feasible options for IRC safe jet 
definitions defined either using sequential recombination [39, 40, 41] or based on the idea 
of cone jets [38] . The recent fast progress in the field of jet physics are expertly reviewed 
in Ref. [42] to which we point the interested reader for further details. 

As a by product of the rapid developments in jet physics there has also recently 
been tremendous interest in using a somewhat more sophisticated understanding of jets 
and their properties, gained via relatively simple analytical calculations, as a chisel for 
improving the prominence of new physics signals at the LHC. For example the idea of the 
optimal value of jet radius R to be employed in various searches for new physics at the 
LHC based on analytical estimates of both perturbative radiation and non-perturbative 
hadronisation corrections was suggested in Ref. [43]. 

Moreover ideas about jet substructure [44, 45] have contributed to an explosion in the 
production of tools which can be used to distinguish QCD jets from those produced by 
the decays of massive particles in the highly boosted regime where the decay products 
may be captured within a single jet. For a detailed exposition of substructure techniques 
we refer the reader to Ref. [46] and references therein. 

To conclude we finish with a reminder that much of the progress in developing QCD 
precision tools and the consequent improvement in understanding QCD effects whether 
in the context of hadron colliders or e+e^ machines should ultimately yield benefits 
beyond the particular context within which it was initiated. For instance the need for 
developing theoretical methods to further the precision that can be obtained via pertur- 
bative techniques at the LHC should in many cases ultimately have spin-offs that would 
pay dividend in the attainment of even higher precision at future linear colliders. There 
is thus much reason to be optimistic in light of the fact that the pace of developments 
of QCD tools continues to be rapid (and possibly even accelerated) stimulated in large 
part by the advent of LHC data. 
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